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Abstract. Around low- and intermediate-mass (1.5-3 M0) red giants, no planets have been found 
inside 0.6 AU. Such a paucity is not seen in the case of IM© main sequence stars. In this study, we 
examine the possibility that short-period planets were engulfed by their host star evolving off the 
main sequence. To do so, we have simulated the orbital evolution of planets, including the effects 
of stellar tide and mass loss, to determine the critical semimajor axis, cicrit, beyond which planets 
survive the RGB expansion of their host star. We have found that flcrit changes drastically around 
2 Mq-. In the lower-mass range, flciit is more than 1 AU, while 

^crit 

small as about 0.2 AU in 

the higher-mass range. Comparison with measured semimajor axes of known planets suggests that 
there is a lack of planets that only planet engulfment never accounts for in the higher-mass range. 
Whether the lack is real affects our understanding of planet formation. Therefore, increasing the 
number of planet samples around evolved intermediate-mass stars is quite meaningful to confirm 
robustness of the lack of planets. 
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INTRODUCTION 

According to recent radial velocity surveys for GK giants, it appears that there is a lack 
of giant planets inside 0.6 AU [18, 13]. Such a paucity is not seen around Sun-like stars 
which often harbor short-period planets such as hot Jupiters. One possibility is that such 
a property is primordial; namely, short-period giant planets (SPGPs) are not present 
around BA dwarfs originally. Another possibility is that SPGPs were removed during 
the evolution of their host star, which is to be explored in this paper. 

A star evolves off the main sequence (MS) towards the red-giant branch (RGB), after 
exhaustion of hydrogen at its center. Then, once the central helium ignites, the star enters 
a next stable phase called the horizontal branch (HB). Most of the giants with detected 
planets are thought to be on their HB. Before reaching the HB, namely, in the RGB 
phase, the star expands substantially, becomes highly luminous, and loses substantial 
mass, which should affect orbits of surrounding planets. 

Villaver & Livio [23], following Sato et al. [18], examined what happened around 
evolving intermediate-mass stars to understand the observed lack of SPGPs around GK 
giants. They demonstrated that Jovian-mass planets initially orbiting at <~ 1 AU around 
2-3Mq stars undergo orbital decay due to stellar tide, ending up being swallowed by their 
host star. However, they made no quantitative comparison with observation. 



The purpose of this paper is to examine whether planet engulfment by host stars is 
responsible for the lack of SPGPs, by making a quantitative comparison between the 
limits of existence of SPGPs predicted theoretically and suggested observationally. To 
this end, we derive a critical semimajor axis beyond which planets survive the host star's 
RGB phase, and investigate its sensitivity to stellar parameters. 

PHYSICAL MODEL 

We simulate the evolution of circular orbit of a planet during its host star's evolution. 
The effects of stellar tide and mass loss on planetary orbit are included. We neglect other 
competing processes such as the frictional and gravitational drag forces by stellar wind 
and change in the planet mass due to stellar- wind accretion and evaporation, which were 
evaluated to be negligible by Villaver & Livio [23] and Duncan & Lissauer [4]. Thus, 
we integrate 

adt TM^\ M^J \ a J M^' ^ ^ 

where a is the semimajor axis, t is time, Mp is the planet's mass. Mi, and 7?* are the mass 
and radius of the host star, respectively, and k is the apsidal motion constant, and T is 
the eddy turnover timescale. The first term on the right-hand side represents the effect 
of stellar tide [7]. Since an RGB star is a slow rotator, we assume no stellar rotation, 
which means that the stellar tide always causes orbital decay of the planet. As for the 
parameters for stellar tide such as k and T, we have followed Rasio et al. [15] and 
Villaver & Livio [23]. 

The second term on the right-hand side of equation (1) represents orbital migration 
due to stellar mass loss (Mi, < 0). We use the Reimers' parameterisation for stellar mass 
loss [16], namely, = — 4 x 10^^^77L^i?^/M^, where r] is the mass-loss parameter; 
is in Mq yr~^ and M^, L^, and i?* are in solar units. 

We simulate stellar evolution directly with the code MESA v2258 [14] to calculate 
Mi,, Ri,, Mi,, and parameters relevant to stellar tide as a function of time t. In this 
calculation in this paper, the effect of convective overshooting is not included (see 
Kunitomo et al. [10] for the results with overshooting). 

ORBITAL EVOLUTION AND SURVIVAL LEVIIT 

Figure 1 shows the orbital evolution of planets around stars with masses of 2Mq (upper 
panel) and 3Mq (lower panel); the metallicity is solar (Z^ = 0.02) in the calculations. The 
solid lines represent the evolution of the stellar radius Ri,. The orbit of a IMj planet with 
different initial semimajor axes (dotted lines) has been integrated from the host-star's 
zero-age main sequence until the planet's semimajor axis falls below the stellar radius 
(a < Ri,) or until the beginning of Asymptotic Giant Branch (AGB) thermal pulse. 

Of special interest in this study is whether planets survive the RGB and HB phases of 
their host star. As seen in the upper panel of Fig. 1, the 2Mq star expands up to 0.5 AU 
at the RGB-tip and consequently swallows planets whose initial semimajor axes, a,, are 
smaller than 1.1 AU. The planet that starts out at 1.2 AU is pulled by the host star, but 
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FIGURE 1. Evolution of semimajor axes of Jupiter-mass planets (dotted lines) and the radius of the 
solar-metallicity star (solid line) for = 2 Mq (upper panel) and 3 Mq (bottom panel). 



it barely survives the RGB/HB phases. We define a critical initial semimajor axis, aan, 
below which planets end up being engulfed by their host star at some point on the RGB 
or HB. In the present case, acrit =1-2 AU. (Hereafter, we sometimes call the critical 
initial semimajor axis the survival limit.) As for the 3Mq star (the lower panel of Fig. 1), 
the RGB-tip radius is as small as 0.13 AU. As a consequence, even a planet starting 
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FIGURE 2. Comparison between the critical semimajor axis acrit (solid line) and measured semimajor 
axes of detected exoplanets (symbols with error bars). The theoretical curve is for stellar metallicity 
Z = 0.01 and planetary mass Mp = 6Mj. In this calculation, the effect of overshooting is not included. 



out at 0.20 AU can avoid engulfment during the RGB and HB phases. As seen in these 
panels, the survival limit is sensitive to the mass of host stars. While acrit = 1-2 AU for 
the 2Mq star, acdt = 0.20 AU for the 3Mq star. 

COMPARISON WITH OBSERVATION 

We compare measured semimajor axes of known planets around GK giants with the 
critical semimajor axis that we have derived in Fig. 2. Those measured stellar masses and 
planetary semimajor axes are taken from the literature [e.g., 22]. Two facts can be found 
in this figure: First, the observational errors being taken into account, it would be fair to 
say that all the known planets are outside the survival limit. Therefore, it is not denied 
that SPGPs (a < acrit) were swallowed by their host star on the RGB. Second, around 
giants of M^^ > 2AMq, planets exist far from acrit- In other words, planet engulfment by 
host stars seems not to be the main reason for the lack of SPGPs. 

Another possibility is that the lack of SPGPs is primordial — there may be any stellar- 
mass -dependent processes that hinder formation of SPGPs around 2-3Mq stars. For 
example, Burkert & Ida [1] proposed that a decreasing-with-stellar-mass lifetime of 
protoplanetary disks could become comparable with or shorter than the timescale of 
the type-II migration of giant planets. They then demonstrated that the observed period 
valley between the hot- Jupiter and cool- Jupiter classes is more pronounced for planets 
orbiting F stars (1.2-1.5M0) than GK stars (O.8-I.2M0). Applying the idea to high-mass 
stars (I.5-3.OM0), Currie [2] demonstrated that SPGPs are rarely formed around high- 



mass stars. 

Obviously we need more samples of giant planets and smaller planets orbiting high- 
mass stars to confirm conclusions we have derived in this paper and to understand stellar- 
mass-dependence of planet formation, which is also helpful in understanding the origins 
and diversity of planetary systems around solar-type stars. Therefore, further surveys for 
planets around giants are highly encouraged. 
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